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Abstract 

A logical error in the usual derivation of the energy conservation law is analyzed, and a 
way to avoid the error is presented. 

In earlier papers [1], [2] we identified a logical error adopted by repetition in textbooks on 
classical electrodynamics when the laws of conservation of energy and momentum are derived 
for a system consisting of electromagnetic fields and charged particles. In this paper we analyze 
the derivation of the energy conservation law, discuss the origin of the logical error made in this 
derivation, and present a way to avoid the error. 

We start by reviewing briefly the usual derivation of the energy conservation law in classical 
electrodynamics. By combining the Maxwell equations in the usual way and integrating over 
a volume V bounded by the surface A, using the divergence theorem, we obtain Poynting's 
theorem [3], [4] 

j t J wdV + I S • hdA + J J • EdV = (1) 

V A V 

In this expression, the quantity 

w = £ i E2 + 2^ 2 ( 2 ) 

is the energy density of the electromagnetic field (SI units are used throughout), E the electric 
field, B the magnetic induction, J = pv the current density, p the charge density, v the velocity 
of the charge in volume element dV, S = E x B//io the Pointing vector, n the unit vector normal 
to the surface A, and dA the element of area of the bounding surface. 

Up to this point there are no inconsistencies, at least for continuous distributions of charge 
and electromagnetic fields that are free of singularities. The problems arise when Poynting's 
theorem is generalized to include point charges, for in this case the total energy in the fields (the 
first term in (1)) diverges and the field E at the position of the charge (in the last term in (1)) 
is not defined. The problems become worse when we generalize (1) to construct a global law for 
the conservation of electromagnetic and mechanical energy, for in this case the self-forces of the 
electromagnetic fields of a particle on itself contribute to the inertia of the particle and to the 
radiation reaction on the particle [3], [4]. We ordinarily include the inertial effect of the self- 
fields in the observed mass of the particle, and therefore double-count it when we add the kinetic 
energy of the particles to the energy of the fields to find the total energy in the system. The 
radiation reaction has its problems as well: it causes unphysical motions of the particle. Thus, 
the electrodynamics of point charges is fraught with contradictions. They persist in quantum 
electrodynamics . 

To simplify the following discussion, we extend the volume V to include all space and consider 
a closed system of particles, so there are no "external" fields coming in from infinity. External 
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fields can be regarded as the fields from other particles that we include in the closed system. 
If the fields vanish sufficiently rapidly at infinity, the integral of the poynting vector over the 
surface A vanishes and we are left with 



dW 
~dT 

where 



+ J J • EdV = 0, (3) 



W = J wdV. (4) 



When the current density arises from a set of point charges qi at positions rj, the current 
density J may be expressed in the form 

J = ^<?iVj<5(r-rj), (5) 

i 

where Vj = dri/dt is the velocity of the i particle. When we substitute this into (3) and 
integrate over all space, the conservation law becomes 

^ + J>E fo) ■ v* = 0. (6) 

i 

But the rate at which work is done on the particle is just the rate of increase of the 
energy ej of the i particle, 

^=«E(r i )-v i . (7) 
If we substitute this into (6), we get the global conservation law 

±(W + E) = 0, (8) 

where 

£ = E e - ( 9 ) 

i 

is the total energy of the particles. 

However, for point charges, this derivation has three problems. In the first place, the electro- 
magnetic energy density w in the fields of the particles diverges at the positions of the particles. 
In the second place, the field E (rj) includes the self-fields Ej (rj). It is not defined at the position 
of the i th particle. Typically, E (rj) is restricted to the field of the other particles, which is well 
behaved at rj, and in the derivation of the energy conservation law, the self- fields Ej (rj) of the 
particle (the inertial and the radiation reaction) are ignored [3], [4]. This is an annoying error. 
In the third place, the energy of the Coulomb field surrounding a particle is generally included 
as part of the mechanical energy of the particle. That is, in the nonrelativistic limit 

ej = ^rriivf, (10) 

where the observed mass mi includes both the "bare mass" of the particle and the kinetic 
energy attributable to the self-electromagnetic field of the particle. Since this electromagnetic 
contribution is infinite, for a point charge, the bare mass of the particle is assumed to be negative 
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and infinite in a way that nearly cancels out the electromagnetic contribution, leaving a finite 
observed mass. This is called renormalization. Even if we set aside questions about the validity 
of dealing with divergent quantities in this way, it still means that we have counted the energy 
of the electromagnetic fields of the particles twice in the conservation law (8), once in W and 
once in E. Thus, we have introduced a logical inconsistency, and the conservation law expressed 
by (7) - (9) is incorrect as it stands. Unfortunately, this error appears in almost every text on 
electrodynamics [3], [4]. Conservation laws for the linear momentum and angular momentum 
(or the 4-vector momentum) similar to (7) - (9) can also be derived, but they have the same 
fundamental problems. Moreover, the energy and momentum of the self-field do not have the 
correct relativistic transformation properties. This is the so-called 4/3 problem, and it is resolved 
only by introducing the so-called Poincare stresses [5], [6], [7] In the following we consider the 
nonrelativistic case, since all the difficulties (including the 4/3 problem) are exhibited, and the 
relativistic effects are not essential to the argument. 

We begin with the last term in (6) or (7), and focus our attention on the self- field contri- 
bution Ej (r-j) to E(rj). To see how to handle the self-interaction, we consider the motion of 
an extended charge distribution and then examine the limit when the size of the distribution 
vanishes. Unfortunately, if we take the limit in the usual way we obtain an equation of motion 
that has nonphysical solutions. The details are described in the Appendix, but for now it is 
enough to point out that the self-electromagnetic force on a small, spherically symmetric, shell 
of charge can be represented by the series 

^(electromagnetic) _ -p , s Qj dvj qf d?Vj 

in the nonrelativistic case, where a is the radius of the charge distribution and brackets \l denote 
averaging through the particle volume. In the nonrelativistic limit, which is valid in the particle 
rest frame, we can ignore magnetic effects. If we overlook the fact that the derivation is based 
on the concept of a rigid charge distribution, which is impossible in relativistic kinematics, 
the relativistic generalization of (11) is straightforward. The first term in (11) is negative 
and proportional to the acceleration. It resists the acceleration in just the way that the bare 
mass does, so we can add it to the inertial term in the equation of motion. It contributes an 
electromagnetic mass 

2 

(electromagnetic) Q 

171 (dynamic) ~ g^~^V ^ > 

which diverges as the particle shrinks to a point (a — > 0). Note that the rest mass that appears 

(electromaqnetic) 4 (electromagnetic) (electromaqnetic) 9 9 /„ , -, r 

is m), \ = ,m; .. -r , where m) .. -r c = q /Sirena is trie energy of 

(dynamic) 6 (energetic) ' (energetic) 1 / u oj 

the electrostatic field around a shell of charge of radius a. This is called the 4/3 problem, and 
the difference is attributable to the so-called Poincare stresses that hold the charged particle 
together [4], [5]. The next term in the expansion (11) of the self- force is called the radiation 
reaction. It remains finite as the particle shrinks to a point, but is responsible for runaway 
solutions, as described in the appendix. We ignore these problems in the following discussion, 
and in fact they can be avoided by taking the limit in a way that represents the self-force by 
a difference-differential equation or integro-differential equation. The remaining terms in the 
series (11) vanish in the limit as the particle shrinks to a point. 

By analogy we can represent the self-nonelectromagnetic force on the same shell in the form 

-p(nonelectromagnetic) (nonelectromagnetic) dv% (13) 

i,(self) ^ i , (dynamic) (it ' 
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The equation (13) is limited to one term (we took into account that the fields of nonelectro- 
magnetic origin are not emitted and hence do not produce a radiation reaction). 

If we postulate that the external force applied to the i th particle is equal to the sum of 
electromagnetic and nonelectromagnetic self-forces and is of opposite direction, then, according 
to (11), (13), the equation of motion can be presented in the form accepted in the classical 
electrodynamics 



m ^served)dv± = pother) + gj d ^ + ( a ), (14) 

at Qirsoc* at z 

-, (observed) (electromaqnetic) , (nonelectromaqnetic) -,11 1 r\(nthpr\ r \ 

where m\ = ra\ {dynami ^ + m\ (dynamic) is the observed mass, E^ ther > (n) = 

(rj) the total electric field at the point rj due to all the other particles j / i, and Ej (r$) 
the field of particle j at the point rj. 

The equation (6) includes the self-fields of the particle Ej(rj). According to (11) and (14) 
the value 

n-V (rAv — n- <f Tf ■ ( r ■ ) ^ v ■ -I- n ■ T?{other) ( r .\ v . _ \ rn {observed) {electromagnetic) ] dVi 

q t iL [r t ) v t -q t < (r t ) > v, + g t i/ (r t ) v, - [m- - «\ ((%riamic) J 

-^{observed) -^{electromagnetic) (T5) 

i i, {dynamic) \ ' 

where j^° hserved ) _ i m {° bserved ) < y 2 ^ - g ^ " ODservec |" kinetic energy of the particles and 

jAelectromaqnetic) 1 {electromaqnetic) 9 , ■, •, . , . •. , ,. , . f , , 

K l{dynamie) = 2 m i,{dyna m ic) v i the kmetlc electromagnetic dynamic energy of the par- 
ticles. 

If we substitute (15) into (6), we obtain the conservation law in the form 

^ \~\y _|_ {({observed) _ ^-{electromagnetic) 

{dynamic) 



^ l (dynamic) 1 

(16) 

where K{° bserved ) = ^{observed) j g ^ e total "observed" kinetic energy of the particles, 

i 

-lAelecromaqnetic) ^ T r{electromaqnetic) ^ -, T r(electromaqnetic) • , , , n , 1 r 

K), . s a ' = yK) ,, . f '. Clearly, K) . f 'is the energy that has been 

(dynamic) t-r 1 i,(dynamic) J ' i,(dynamic) OJ 

double counted in (8) [1],[2] . 

The equation (16) still contains divergent terms in W and ^ e J ecroma 9 net ' LC ) go ^ j s impossible 

^ v / (dynamic) 1 

to use the conservation law in this form. To deal with this, we can proceed in the following 
manner. Since the field near the particle approaches the Coulomb field of a homogeneously 
moving particle, we can write 

-r-, -^{Coulomb) -^{remainder) r-tn\ 

&i = tit , (17) 

with a similar decomposition of the magnetic field. In the limit as the size of the particle 
vanishes, we see that its electromagnetic energy is 

_ ^{Coulomb) _|_ -^{remainder) 

where the reminder energy of the \ th particle is 

(remainder) ^0 J ^^{Coulomb) -^{remainder) ^ -^{remainder)' 2 ^ 

^1 f ( 2g(C 'oulomb) ^{remainder) ^ j^{remainder) 2 \ (19) 

J \ ) 
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which is finite. 

In this case we can represent the total energy of the electromagnetic fields in the form 



W = V. Wi + § / V Ei-EjdV + — [y^Bi- BjdV = w Coulomb + VF re ~ der , (20) 

where the value w( Coulornb ) = £\ ^( Coulomb ) corresponds to the energy of particles homoge- 
neously moving with the given velocity at the moment of observation and the electromagnetic 
energy w { - remainder ^ = wf emainder ^ + f / £ Ei-EjdV +^/EB r BjdV is the sum of the 

interaction energy of particles with electromagnetic fields and the energy of emitted radiation. 
The energy in the self-field of a homogeneously moving spherical shell of a charge is 

w (Coui omb) = so J E jdV + ±- J B^dV = m ^;:~^ 

r>a r>a 



(electromagnetic) 2 , 5 ^(electromagnetic) 2 (electromagnetic) 2 1 5 -^(electromagnetic) (21) 

^i, (energetic) ^ 6 i, (energetic) ^i, (energetic) ^ 3 i, (energetic) ' v / 

1 jy- (electromagnetic) 1 (electromagnetic) 2 / rol\ 

WUere ^■ ii ( enerfletic ) - 2 m i,(energetic) ^ V see > e -S-> l 8 JJ- 

If we substitute the total electromagnetic energy of the system W = w^ remainde ^+W^ Coulom ^ 
to (16) and neglect the derivative of the constant terms ^^neVgeUc)^ ' ^ ^ * ne ec t ua ti° n takes 
the form 



d_ 
~dt 



yyr (remainder) , ^(observed) , {^-^(electromagnetic) _ j^(electromagnetic)\ 

V 3 (energetic) (dynamic) J 



0, (22) 



where ) = EW f<— ) and = K^^Uc)^ is the kinetic 

electromagnetic energy of particles homogeneously moving with the velocity Vj. 

The difference term in parentheses of (22) is ^ner^tic)^^ 1^ ^ 0- We can su PP ose 
that this difference term is attributable to the so-called Poincare stresses that hold the charged 
particle together, that if the electromagnetic energy is transformed by a complicated law (21), the 
nonelectromagnetic energy is transformed by any analogous law and compensate this difference. 
If we ignore it then this leaves the conservation law 



d_ 

dt 



^■(remainder) _|_ ^(observed) J q (23) 

The conservation law (23) can be expressed in the equivalent form 

d -^(nonelecromagnetic) 



dt I" + K (energetic) * 'J -0 ' ( 24 ) 

h (nonelectromagnetic) -^(nonelectromagnetic) \ ^(nonelectromagnetic) -^(nonelectromagnetic) 

(energetic) (dynamic) ' 3 (energetic) ' (dynamic) 

Y~\ 1 (nonelectromagnetic) 2 
t-r< 2 m i, (dynamic) V i " 

In this form, the conservation law (24) states explicitly that the sum of the electromag- 
netic energy and the kinetic energy attributable to the nonelectromagnetic energetic energy of 
particles is a constant. It still contains divergent terms. We can extract from W divergent 
terms corresponding to the accompanying electromagnetic energy of particles, combine them 
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with K^°™f g ^™ magnettc ' > and postulate that the obtained value is the observed energy of parti- 
cles. This can be called renormalization. It leaves the conservation law (23) in the same form 
but now the observed kinetic energy is composed of electromagnetic and nonelectromagnetic 
energies determined by masses of energetic origin. 

Both the electromagnetic and nonelectromagnetic dynamic masses in the equation of motion 
(14) can be considered as coefficients at the resistive terms. They have the dimensions of mass 
and their sum agrees with the observable mass. The sum of electromagnetic and nonelectro- 
magnetic energetic masses of particles is another form of presentation of the observable mass of 
particles based on their energy. 

There is still the matter of the 4/3 problem and the energy due to the Poincare stresses. In 
the nonrelativistic theory we simply ignore the divergent, constant terms in the field energy W, 
as discussed above, so the problem disappears. In the relativistic case, however, this energy is 
part of the rest energy of the particle. Unfortunately, this energy is only 3/4 what we would 
expect from the mass that appears in the momentum. The missing mass, as first pointed out by 
Poincare, is accounted for by considering the forces that hold the electric charge distribution. In 
the simplest example, the momentum density of the Poincare stresses (the off-diagonal elements 
of the stress 4-tensor) vanish in all coordinate systems, and the energy density is just sufficient 
to make up the missing mass. More elaborate models of the Poincare stress tensor have both 
momentum and energy, but when the Poincare stresses are included in the symmetric stress 
tensor, a covariant form of the conservation law can be derived [4], [5]. 

To conclude, we have shown a logical inconsistency in the derivation of the energy conserva- 
tion law that appears in almost every text on advanced electrodynamics. The conservation law 
expressed by (7) - (9) is incorrect as it stands and must be replaced by (23). Unfortunately we 
can not state that the correct formulae for the energy conservation law was obtained without any 
hypotheses. It is better to say that we proceed from the assumption that the energy conservation 
law and Maxwell equations are valid in any case and put limitations on the nonelectromagnetic 
fields. 

Appendix 

We start by reviewing very briefly the derivation of the rate of change of momentum of the 
particle based on Abraham-Lorentz evaluation of the self-force. In the nonrelativistic case, the 
equation of motion of a particle in the external fields can be written in the form 

m {bare)<^_ = g pother) + y x ^(other)^ + F (self) ^ ^ 

where 

m (bare) j g 

the mass of the bare particle (without electromagnetic fields), ~E,( other ) and 
-g(other) are ^ e electric and magnetic fields of other particles in the system, and F^ sel ^ is the 
electromagnetic force of the charge distribution of the particle back on itself. In the nonrela- 
tivistic limit we can ignore the effect of the self-magnetic field compared to that of the electric 
field, and the effect of the self-electric field can be expressed as an integral over the retarded 
self-field of the particle of the form 



F( "" ,(t) - ikj d3 ^> I ^' m ( R a? " iw) + w> 



1 (~d P ldJ\ R 



(26) 

retarded 



where 

R = r — r', (27) 
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and the retarded time is 



|r — r" 



t — t re tarded — t . (28) 

For a rigid charge distribution, the charge density at the retarded time t re t is simply related 
to that at the present time t by the motion of the center of mass of the particle. If the charge 
distribution is small, the retarded times are all close to the present time and we can use a 
Taylor-series expansion to evaluate p(r',t re t). For a rigid, spherical shell of charge the result is 

F W) = ? y (zlT (*Y m 

12vre ca 3 ^ (n)! \cj dt n ' 1 ' 

where q is the total charge, a the radius, and v (t) the velocity of the charge [4], [8]. The series 
can be summed, and the result is given by the expression 



2a - 

— -v(t) 



(30) 



l2ir£oca 2 

When this is substituted into the nonrelativistic equation of motion, we get the difference- 
differential equation 



'— = a (^ other ) + v x B(° ther A + q2 ^- ^ 
dt H \ " ' ~ / 127re n a 2 



c 



(31) 



This is called the Page-Somerfeld equation of motion [9], [10]. Its relativistic generalization is 
discussed by Caldirola [11]. For more general charge distributions we get an integro-differential 
equation of motion called the Markov equation [12]. Its relativistic generalization is discussed 
by Brau [7]. 

If instead of summing the series (29) we take just the first two terms, we get 

F (self) = Q 2 dv + q 2 d 2 v^ 
67reoac 2 dt 67reoc 3 dt 2 
If we substitute this into the equation of motion (25) we get 

m (observed)dv = (pother) + v x B (other)\ + (33) 

dt V / 67reoa dt 2 

where the observed mass of the particle is given by (14). The result (33) is called the Abraham- 
Lorentz equation of motion [13], [14]. Its relativistic generalization is discussed by Dirac [15]. 
Unfortunately, this equation admits runaway solutions, for if the external fields vanish the 
equation of motion (33) is satisfied by the solution 

v(t)=v e*/ T , (34) 

where r = 2a /3c and vq is a constant. That is, in the absence of external fields the particle can 
start at rest and accelerate without limit. Runaway solutions are avoided by the Page-Somerfeld 
equation of motion and, under certain conditions, by the Markov equation of motion. However, 
these equations of motion admit oscillatory solutions [16]. 
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